Isotherm, Thermodynamic and Kinetic Studies of Elemental Sulfur Removal from Mineral Insulating Oils Using Highly Selective Adsorbent

Elemental sulfur (S8) is a corrosive sulfur compound which was found to be extremely reactive to silver, causing intensive silver sulfide (Ag2S) deposition on on-load tap changer (OLTC) contacts in power transformers. A highly selective adsorbent (HSA), called Tesla’Ssorb, for the removal of S8 from mineral insulating oils was prepared from raw material (RM) using the novel procedure. In this study, the adsorption properties of HSA for the removal of S8 from the oil were determined. RM and HSA were characterized using various techniques, such as field-emission scanning electron microscopy (FESEM), energy-dispersive X-ray (EDX), and X-ray diffraction (XRD). The performance of HSA was determined by adsorption equilibrium, thermodynamic, and kinetic study through batch experiments, at various temperatures and initial concentrations of S8. The obtained results were analyzed by Langmuir and Freundlich adsorption isotherms and it was found that equilibrium data were fitted better with the Langmuir isotherm model. The maximum adsorption capacity was 4.84 mg of S8/g of HSA at 353 K. Thermodynamic parameters, such as enthalpy (ΔH°), Gibbs free energy (ΔG°), and entropy (ΔS°), were calculated and it was found that the sorption process was spontaneous (ΔG° < 0) and endothermic in nature (ΔH° > 0). It was found that the adsorption of S8 follows pseudo-second-order kinetic model, and the activation energy indicated the activated chemisorption process.


Introduction
The presence of corrosive sulfur in mineral insulating oils can be the root cause of power transformers and on-load tap changer (OLTC) failures, owing to the formation of electro-conductive sulfides in reaction with power transformer construction metals, copper and silver. The mechanisms of copper sulfide formation from di-benzyl-disulfide (DBDS) have been thoroughly investigated over last decades [1][2][3].
In recent times, cases of copper or silver sulfides formation after oil reclamation have been reported [3][4][5][6][7]. Moreover, the large number of OLTC failures recently reported was connected with the corrosion of silver-plated contacts, after the oil reclamation process, conducted to remove oil ageing products or DBDS. Elemental sulfur (S 8 ) was confirmed to be the cause of failure [3,8,9]. In almost all of the cases reported, when silver corrosion was severe and caused by S 8 , copper surfaces were clean, without copper sulfide deposits [4][5][6]. S 8 may be created as a consequence of the application of used mineral insulating oil regeneration processes with adsorbents, as one of the products of the high-temperature combustion process, during reactivation of the adsorbent. If it remains in oil after reclamation, S 8 will easily react with OLTC's silver coated contacts at lower operating temperatures and form conductive silver sulfide deposits during transformer service [1]. Silver sulfide semi-conductive deposits will increase in contact with resistance which will further lead to overheating the OLTC's contacts. Further, once flaked off and suspended in the transformer, there will be an increase in terms of oil conductivity and the electric stress between the contacts [3,10,11]. Consequently, the dielectric strength of oil will be reduced. Furthermore, electric discharge of the OLTC tap selector contacts in the main tank of the transformer may occur, which may cause a short circuit in the regulating winding [12,13]. The typical concentrations of S 8 , found in mineral oils in service, were up to 10-15 mg/kg [1]. Due to high affinity of S 8 to react with silver and a large amount of oil in the transformer compared to the small surface of OLTC's silver plated contacts, only a few mg/kg of S 8 in the oil is sufficient to create silver corrosion, at low operating temperatures (below 60 • C) [10,14].
One of the most frequently applied temporary mitigation techniques to suppress copper corrosion and the deposition of copper sulfide in the transformer windings is the use of metal passivators, but this has been found to be inefficient for the mitigation of silver corrosion [1, [15][16][17].
The mitigation method which provides a permanent solution for the corrosive sulfur problem is oil treatment. It was observed by experience that S 8 is more complex to remove from the oil than DBDS, which is also corrosive to copper and silver [1,4,5,[18][19][20][21]. The solvent extraction process has been proven to be efficient in the removal of corrosive sulfur compounds with a significant reduction of aromatic compounds as well, indicating improvement of oil properties and degree of oil refining [22][23][24][25]. Desulfurization processes with K-PEG or Na-PEG reagent are efficient in the removal of DBDS from the oil, but not for S 8 removal [26].
Oil treatment processes with adsorbents are based on adsorption or chemisorption of unwanted compounds from mineral oil on the adsorbent surface, by forced circulation of oil thought columns with packed adsorbent bed. Conventional adsorbents, such as aluminosilicates (Fuller's earth), magnesium aluminium silicate (Ultrasorb), and aluminium oxides (Bauxite), are widely used in oil reclamation processes. They are found to be very efficient in the removal of polar oil ageing products and some corrosive sulfur compounds, namely DBDS. In the case of S 8 removal, these adsorbents are found to be mostly completely inneficient in the removal of S 8 , especially when S 8 is present in high concentrations, or of low efficiency when S 8 was present in low concentrations (up to 10 mg/kg of S 8 ) [1,4].
HSA is a tailor-made adsorbent, called Tesla'Ssorb, based on silicon dioxide with bonded silver, specially designed to selectively remove reactive sulfur compound S 8 (in chemical reaction with deposited silver) in low amounts (up to 3 wt.% of adsorbent) and in short treatment time, in comparison to other adsorbents. The previous study indicated the efficiency of HSA to remove S 8, DBDS and oil ageing products, obtaining the oil with restored properties for further use in electrical equipment [27]. Moreover, it was shown that HSA was efficient in the removal of high concentrations of S 8 (approximately 45 mg/kg) to a final concentration below 1 mg/kg of S 8 , and therefore it can be successfully used to reduce the risks of power transformers failures and to extend transformer life [27].
In this study, the adsorption performance of Tesla'Ssorb by treatment the oil at various temperatures and initial concentrations of S 8 was investigated. The adsorption process was studied via isotherm, thermodynamic, and kinetic studies. Characterization of the HSA was performed using various analytical methods, in order to define the morphology, phase, and elemental composition of the HSA and to comprehend the adsorption mechanism better.
The kinetic of the adsorption process was studied at three different temperatures, 328, 338, and 353 K. The mechanism of binding of S 8 molecules to the HSA's surface was revealed by determining the activation energy for the process from the kinetics data using Arrhenius plots.

Materials and Reagents
Tesla'Ssorb was synthesized by the impregnation method using an adsorbent containing pure silicon dioxide (SiO 2 , 94-98 wt.%) and traces of calcium oxide (CaO, 0.04-1 wt.%) as raw material (RM). The chemicals used for the synthesis of HSA were as follows: silver nitrate (AgNO 3 , extra pure >99.98%, Bor, Serbia), aqueous ammonium hydroxide (25%, Macron, Austria), and water, deionized, H 2 O (Pharma Product, Belgrade, Serbia). Sulfur, powder, reagent grade was procured from Sigma Aldrich (St. Louis, MO, USA), while di phenyl disulfide (DPDS, ≥99%) was obtained from Merck, Germany. The standard solutions and formulated oil samples (spiked with known amount of S 8 ) were prepared in used mineral insulating oil free from S 8 and DPDS. The preparation of oil samples was done with toluene and iso-octane. All chemicals used were of analytical grade and without further purification.

Adsorbent Preparation and Characterization
Tesla'Ssorb for the removal of S 8 from mineral insulating oils was synthesized by the impregnation method. The procedure comprised of three stages:

1.
Annealing the adsorbent support consisted mainly of silicon dioxide (RM) at 150 • C, for 18-24 h, to eliminate adsorbed moisture.

2.
Deposition of silver ions (3-6 wt.%) on the adsorbent support by silver nitrate aqueous solution followed by annealing at temperature from 120 • C to 130 • C, during 18-24 h.

3.
Deposition of ammonium hydroxide aqueous solution (5-10 wt.%) on the adsorbent support followed by annealing the adsorbent at temperature from 125 • C to 130 • C, during 18-24 h. Ammonium hydroxide is added to neutralize the acidic by-products formed during removal of S 8 from the oil [28].
For producing an adsorbent with high effectiveness for the removal of sulfur compounds corrosive to silver and to obtain restored oil properties with low acidity for further use in power transformers, both compounds, i.e., silver nitrate and ammonium hydroxide, are necessary [27].
The removal of S 8 is obtained by the chemical reaction with deposited silver ions on the adsorbent surface, through the chemisorption process, followed by the addition of ammonium hydroxide for the neutralization of acidic by-products. This is shown in the following reaction 1: An effective binding of silver ions to Tesla'Ssorb silanol groups was achieved with this temperature induced impregnation method. Moreover, the formation of silver oxides was minimized, which may reduce reactivity of silver with corrosive sulfur compounds in the oil. Furthermore, the addition of ammonium hydroxide in a specific and optimal concentration range is necessary for the efficient neutralization of acidic by-products. On the other side, the said compound is not present in an excess amount, since it can reduce the mobility of silver ions for the reaction and the efficiency of corrosive sulfur removal, due to the possible formation of complex salts with silver, i.e., diamine silver (I) [(Ag(NH 3 ) 2 ] + complex [27].
Basic physical properties of Tesla'Ssorb were analyzed, applying the standard procedures (ASTM, EPA) and in-house validated methods. Field-emission scanning electron microscopy (FESEM) (TESCAN MIRA3XMU, Brno, Czech Republic) operated at 20 keV was used to observe the surface morphology of RM and HSA. An atomic gold layer was deposited on the sample surfaces before analysis. Moreover, energy-dispersive X-ray (EDX) (Oxford Inca 3.2 coupled to a JEOL JSM 5800 scanning electron microscope (JEOL, Akishima, Tokyo, Japan)) analysis was used to analyze the elemental composition of a synthesized adsorbent. X-ray diffraction (Rigaku MiniFlex600 equipped with D/teX Ultra 250 high-speed detector) was performed using CuKα (λ = 1.54 Å), operated at 40 kV and 15 mA, in the 2θ angular range of 3-90 • , with 0.02 • of step size and 10 • /min of scan speed. Measurements were performed using the software MiniFlex Guidance Version 2.1.0.4. The mineral identification was performed in the software PDXL 2 Version 2.4.2.0 and the obtained diffractograms were compared with the data from the database ICDD (PDF-2 Release 2015 RDB).

Preparation and Analysis of S 8 in Mineral Insulating Oils
Solutions of S 8 were prepared by dissolving a necessary amount of sulfur powder in used mineral insulating oil free from S 8 , in order to obtain the desired concentrations of S 8 ranging from approximately 80 to 150 mg/kg. The concentrations of S 8 in used mineral oils were determined using GC-ECD (Thermo Scientific 1300 system with an autosampler AI 3000 series, Milan, Italy), according to IEC TR 62697-3.

Batch Adsorption Experiments
The batch adsorption experiments were evaluated in terms of the effect of initial S 8 concentrations, contact time, and temperature, in order to study adsorption isotherms, adsorption kinetics, and thermodynamic. Experiments were carried out using the percolation process, on a pilot scale reclamation unit with a digital temperature controlled hot plate ( Figure 1). The 17.3 kg of used mineral insulating oil, spiked with various initial concentrations of S 8 from approximately 80 to 150 mg/kg was pumped through a stainless-steel column filled with 523 gr of HSA (i.e., 3 wt.% of adsorbents compared to the mass of oil) at a flow rate of ca. 200 mL/min. To achieve equilibrium at constant temperatures of 328, 338, and 353 K respectively, a contact time of 1000 min was fixed in all experiments. At predetermined times, the samples were collected and analyzed. microscopy (FESEM) (TESCAN MIRA3XMU, Brno, Czech Republic) operated at 20 keV was used to observe the surface morphology of RM and HSA. An atomic gold layer was deposited on the sample surfaces before analysis. Moreover, energy-dispersive X-ray (EDX) (Oxford Inca 3.2 coupled to a JEOL JSM 5800 scanning electron microscope (JEOL, Akishima, Tokyo, Japan)) analysis was used to analyze the elemental composition of a synthesized adsorbent. X-ray diffraction (Rigaku MiniFlex600 equipped with D/teX Ultra 250 high-speed detector) was performed using CuKα (λ = 1.54 Å), operated at 40 kV and 15 mA, in the 2θ angular range of 3-90°, with 0.02° of step size and 10°/min of scan speed. Measurements were performed using the software MiniFlex Guidance Version 2.1.0.4. The mineral identification was performed in the software PDXL 2 Version 2.4.2.0 and the obtained diffractograms were compared with the data from the database ICDD (PDF-2 Release 2015 RDB).

Preparation and Analysis of S8 in Mineral Insulating Oils
Solutions of S8 were prepared by dissolving a necessary amount of sulfur powder in used mineral insulating oil free from S8, in order to obtain the desired concentrations of S8 ranging from approximately 80 to 150 mg/kg. The concentrations of S8 in used mineral oils were determined using GC-ECD (Thermo Scientific 1300 system with an autosampler AI 3000 series, Milan, Italy), according to IEC TR 62697-3.

Batch Adsorption Experiments
The batch adsorption experiments were evaluated in terms of the effect of initial S8 concentrations, contact time, and temperature, in order to study adsorption isotherms, adsorption kinetics, and thermodynamic. Experiments were carried out using the percolation process, on a pilot scale reclamation unit with a digital temperature controlled hot plate ( Figure 1). The 17.3 kg of used mineral insulating oil, spiked with various initial concentrations of S8 from approximately 80 to 150 mg/kg was pumped through a stainlesssteel column filled with 523 gr of HSA (i.e., 3 wt.% of adsorbents compared to the mass of oil) at a flow rate of ca. 200 mL/min. To achieve equilibrium at constant temperatures of 328, 338, and 353 K respectively, a contact time of 1000 min was fixed in all experiments. At predetermined times, the samples were collected and analyzed.  The percentage of S 8 adsorption (%) and adsorption capacities (q e and q t ) were calculated using Equations (2)-(4): where q e is the amount of adsorbate adsorbed at equilibrium per unit weight of the adsorbent (mg g −1 ), q t is the amount of adsorbate adsorbed at any time (mg g −1 ), C 0 and C e are the initial and equilibrium concentrations of S 8 (mg·kg −1 ), m oil is the mass of oil (kg), and m ads is the adsorbent mass (g). In order to determine the kinetic model, limiting level, and mechanism of the adsorption of S 8 onto the HSA, the pseudo-first-order [29], pseudo-second-order [30], and intraparticle diffusion kinetic models were used [31]. The Langmuir [32] and Freundlich [33] models were used to describe the equilibrium properties of S 8 adsorption on HSA. The thermodynamic parameters, such as the change of Gibbs free energy (∆G o • ), enthalpy (∆H o • ), and entropy (∆S o • ) were calculated in order to confirm the reaction mechanism (reaction 1). Insight into the binding process was provided by the determination of activation energy using the Arrhenius plot.

Characterization of Tesla'Ssorb
The Tesla'Ssorb was characterized by determining some of its general physical properties, e.g., pH, flash point, melting point, and bulk density, employing standard methods. The results are shown in Table 1. Based on the results of textural characteristics of HSA, such as specific surface area, S BET (29.28 m 2 /g), total pore volume, V total (0.2026 cm 3 /g), mesopore volume, V meso (0.2025 cm 3 /g), average pore diameter, D A (19.2 nm), and pore size distribution, D BHJ (25.1 nm), reported in previous work, HSA was classified as a mesoporous material [27].
The FESEM micrographs of RM (before) and Tesla'Ssorb (after) modification, i.e., impregnation and sorption of Ag ions, are displayed in Figures 2 and 3.
Based on the FESEM surface analysis of the RM (before) and Tesla'Ssorb (after) sorption of Ag ions at lower magnifications, irregularly shaped agglomerates with dimensions larger than 1 mm are observed (Figure 2a). At higher magnifications, it was observed that both RM and Tesla'Ssorb agglomerates consist of rod-like particles and particles of irregular shape (Figures 2b,c and 3b). In the case of the Tesla'Ssorb, rod-like particles were found to be more elongated than the primary RM particles. Further analysis of Tesla'Ssorb at the highest magnifications shows the existence of nano-sized spherical clusters on the surface of the primary particles ( Figure 3c). This phenomenon has not been observed in the morphology of the RM (Figure 2c). Therefore, it can be concluded that these spherical clusters, with a diameter of approximately 60 nm, are silver nanoparticles formed during the modification of the RM, by silver nitrate impregnation. Materials 2023, 16, x FOR PEER REVIEW 6 of 21  Based on the FESEM surface analysis of the RM (before) and Tesla'Ssorb (after) sorption of Ag ions at lower magnifications, irregularly shaped agglomerates with dimensions larger than 1 mm are observed (Figure 2a). At higher magnifications, it was observed that both RM and Tesla'Ssorb agglomerates consist of rod-like particles and particles of irregular shape (Figures 2b,c and 3b). In the case of the Tesla'Ssorb, rod-like particles were found to be more elongated than the primary RM particles. Further analysis of Tesla'Ssorb at the highest magnifications shows the existence of nano-sized spherical clusters on the surface of the primary particles ( Figure 3c). This phenomenon has not been observed in the morphology of the RM (Figure 2c). Therefore, it can be concluded that these spherical clusters, with a diameter of approximately 60 nm, are silver nanoparticles formed during the modification of the RM, by silver nitrate impregnation.
The X-ray diffraction patterns of RM and Tesla'Ssorb are shown in Figure 4a   Based on the FESEM surface analysis of the RM (before) and Tesla'Ssorb (after) sorption of Ag ions at lower magnifications, irregularly shaped agglomerates with dimensions larger than 1 mm are observed (Figure 2a). At higher magnifications, it was observed that both RM and Tesla'Ssorb agglomerates consist of rod-like particles and particles of irregular shape (Figures 2b,c and 3b). In the case of the Tesla'Ssorb, rod-like particles were found to be more elongated than the primary RM particles. Further analysis of Tesla'Ssorb at the highest magnifications shows the existence of nano-sized spherical clusters on the surface of the primary particles (Figure 3c). This phenomenon has not been observed in the morphology of the RM (Figure 2c). Therefore, it can be concluded that these spherical clusters, with a diameter of approximately 60 nm, are silver nanoparticles formed during the modification of the RM, by silver nitrate impregnation.

Effect of Initial S 8 Concentration
The effect of initial S 8 concentration in the range of 79.7-153.1 mg/kg on adsorption using Tesla'Ssorb was carried out at different temperatures of 328, 338, and 353 K, and the results are given in Figure 6.

Effect of Initial S8 Concentration
The effect of initial S8 concentration in the range of 79.7-153.1 mg/kg on adsorption using Tesla'Ssorb was carried out at different temperatures of 328, 338, and 353 K, and the results are given in Figure 6. Due to the fixed amount of adsorbent used in this study (3 wt.%), it can be seen from the Figure 6 that the percentage of S8 adsorption decreased as the initial concentration of S8 increased. The presence of more molecules of S8 in the oil per unit number of adsorbent sites (with bonded silver ions) leads to a saturation of the adsorption sites.
For the 79.7 mg/kg of S8, the percentage of adsorption was found to be approximately 99.0% at all temperatures, while for 153.1 mg/kg of S8, the relevant values were 83.2%, 84.8%, and 91.4% at 328, 338, and 353 K respectively.
Moreover, the higher amount of S8 is adsorbed by Tesla'Ssorb at the highest initial concentrations of S8 (153.1 mg/kg). The initial concentration of S8 has a significant impact on the adsorption capacity since the concentration gradient provides the necessary driving force for adsorption between the oil and the Tesla'Ssorb surface.
The increase in initial S8 concentration improved the interaction between incorporated silver ions on the Tesla'Ssorb and S8 molecules as well as the adsorption capability of Tesla'Ssorb to remove S8 from the oil, even if S8 is present in the oil in very high concentrations. Figure 7 illustrates the relation between contact time and S8 adsorption on Tesla'Ssorb at various initial S8 concentrations at highest temperature (353 K). The adsorption was very fast during the first 300 min. With a further increase of time, the rate of adsorption decreased, and the adsorption equilibrium was achieved within 600 min, i.e., within 1000 min (for the highest concentration of S8, 153.1 mg/kg), as shown in Figure 7. Due to the fixed amount of adsorbent used in this study (3 wt.%), it can be seen from the Figure 6 that the percentage of S 8 adsorption decreased as the initial concentration of S 8 increased. The presence of more molecules of S 8 in the oil per unit number of adsorbent sites (with bonded silver ions) leads to a saturation of the adsorption sites.

Effect of Contact Time
For the 79.7 mg/kg of S 8 , the percentage of adsorption was found to be approximately 99.0% at all temperatures, while for 153.1 mg/kg of S 8 , the relevant values were 83.2%, 84.8%, and 91.4% at 328, 338, and 353 K respectively.
Moreover, the higher amount of S 8 is adsorbed by Tesla'Ssorb at the highest initial concentrations of S 8 (153.1 mg/kg). The initial concentration of S 8 has a significant impact on the adsorption capacity since the concentration gradient provides the necessary driving force for adsorption between the oil and the Tesla'Ssorb surface.
The increase in initial S 8 concentration improved the interaction between incorporated silver ions on the Tesla'Ssorb and S 8 molecules as well as the adsorption capability of Tesla'Ssorb to remove S 8 from the oil, even if S 8 is present in the oil in very high concentrations. At adsorption equilibrium, a constant value of S8 is obtained, where no more S8 is removed from the oil. These results are in very good correlation with the equilibrium adsorption capacity of Tesla'Ssorb previously reported, for the treatment of oil from a real 35 kV power transformer, performed on-site [27].

Effect of Contact Time
The fast initial adsorption is due to the increased concentration gradient between the S8 in the oil and incorporated silver ions on the adsorbent support, since there are a large number of unoccupied sites available initially. Achieving the adsorption equilibrium over time is due to the reduced number of available adsorbent sites for reaction with S8 molecules and the low concentration of S8 in the oil. At the highest adsorption temperature of 353 K, the percentage of adsorption at equilibrium stage decreased from 99.1% to 91.4% with increasing S8 concentration from 79.7 to 153.1 mg/kg.

Effect of Temperature
The removal of S8 using Tesla'Ssorb was studied as a function of temperature, with experiments on oils containing different initial S8 concentrations of (79.7 to 153.1) mg/kg, at different temperatures of (328, 338, and 353) K using a fixed adsorbent dosage of 3 wt.% relative to the mass of oil (Figure 8).
At the lowest temperature (328 K), a significant difference in percentage of adsorption was observed, for various initial concentrations of S8 (79.7, 135.7 and 153.1) mg/kg in comparison with results obtained at 353 K, where the degree of S8 removal was more than 90% (from 91.4% to 99.1%), for the same initial concentrations of S8.
These results indicated that higher temperatures increase the activity of incorporated silver ions on the Tesla'Ssorb to react with S8 molecules present in the oil, promoting the removal of S8 from the oil in a wide range of concentrations. Therefore, it can be concluded that 353 K can be used as an optimal temperature for S8 removal for the oil, in a wide range of concentrations. At adsorption equilibrium, a constant value of S 8 is obtained, where no more S 8 is removed from the oil. These results are in very good correlation with the equilibrium adsorption capacity of Tesla'Ssorb previously reported, for the treatment of oil from a real 35 kV power transformer, performed on-site [27].
The fast initial adsorption is due to the increased concentration gradient between the S 8 in the oil and incorporated silver ions on the adsorbent support, since there are a large number of unoccupied sites available initially. Achieving the adsorption equilibrium over time is due to the reduced number of available adsorbent sites for reaction with S 8 molecules and the low concentration of S 8 in the oil. At the highest adsorption temperature of 353 K, the percentage of adsorption at equilibrium stage decreased from 99.1% to 91.4% with increasing S 8 concentration from 79.7 to 153.1 mg/kg.

Effect of Temperature
The removal of S 8 using Tesla'Ssorb was studied as a function of temperature, with experiments on oils containing different initial S 8 concentrations of (79.7 to 153.1) mg/kg, at different temperatures of (328, 338, and 353) K using a fixed adsorbent dosage of 3 wt.% relative to the mass of oil ( Figure 8). Results presented herein showed that the adsorption capacity of Tesla'Ssorb was enhanced as the temperature was increased, indicating the endothermic nature of the S8 adsorption process. At the lowest temperature (328 K), a significant difference in percentage of adsorption was observed, for various initial concentrations of S 8 (79.7, 135.7 and 153.1) mg/kg in comparison with results obtained at 353 K, where the degree of S 8 removal was more than 90% (from 91.4% to 99.1%), for the same initial concentrations of S 8 .

Adsorption Isotherm Study
These results indicated that higher temperatures increase the activity of incorporated silver ions on the Tesla'Ssorb to react with S 8 molecules present in the oil, promoting the removal of S 8 from the oil in a wide range of concentrations. Therefore, it can be concluded that 353 K can be used as an optimal temperature for S 8 removal for the oil, in a wide range of concentrations.
Results presented herein showed that the adsorption capacity of Tesla'Ssorb was enhanced as the temperature was increased, indicating the endothermic nature of the S 8 adsorption process.

Adsorption Isotherm Study
The experimental data were analyzed using Langmuir and Freundlich isotherms to fully understand the type of interaction between the adsorbate and the adsorbent. The Langmuir isotherm model is based on the premise that the adsorbent is structurally homogeneous, covered in a monolayer with no interaction between the molecules of the absorbate. The equation developed by Langmuir is expressed as follows: where q e is the amount of S 8 adsorbed at equilibrium per unit weight of the adsorbent (mg·g −1 ), C e is the equilibrium concentrations of S 8 (mg·kg −1 ) in the oil, K L is the Langmuir constant (kg·mg −1 ), and q m is the maximum adsorption capacity (mg·g −1 ). The maximum adsorption capacity (q m ) of Tesla'Ssorb and Langmuir constant were determined from the curves C e /q e vs. 1/C e (Figures 9-11). It was found that the maximum adsorption capacity of Tesla'Ssorb is 4.84 mg·g −1 , obtained at 353 K ( Table 2). The separation factor (RL) can be used to describe the Langmuir isotherm's feasibility and may be calculated using the following equation: The RL value defines the form of isotherm and the adsorption process nature, which is irreversible (RL = 0), unfavorable (RL > 1), linear (RL = 1), and favorable (0 < RL < 1). The obtained RL values in this study were in the range of 0−1, thus confirming that the Tesla'Ssorb is a suitable adsorbent for the adsorption of S8 under defined testing conditions.
Graphical presentations for Langmuir and Freundlich equilibrium isotherms of S8 adsorption onto Tesla'Ssorb (linear and non-linear models) are given in Figures 9-11. Correlation coefficients (R 2 values) were used to compare the suitability of used isotherm models. Since the obtained R 2 values for Langmuir isotherm model were higher in comparison with Freundlich isotherm (Table 2), it can be concluded that the Langmuir isotherm fits better. The experimental data, in both linearized and non-linearized fitting curves, were fitted better to the Langmuir isotherm model, as shown in Figures 9-11. Furthermore, it can be observed that maximum adsorption capacity (qm) increases with increasing temperature, due to the endothermic nature of adsorption. This was further confirmed by the analysis of thermodynamic parameters.     In comparison with Langmuir isotherm model, the Freundlich isotherm presumes the interaction between adsorbed molecules on heterogeneous surfaces. This isotherm is expressed as follows: where K F and n are Freundlich constants related to the adsorption capacity and adsorption intensity, respectively. The typical parameters for the Langmuir and Freundlich isotherm models and R 2 values at three temperatures are given in Table 2.
The separation factor (R L ) can be used to describe the Langmuir isotherm's feasibility and may be calculated using the following equation: The R L value defines the form of isotherm and the adsorption process nature, which is irreversible (R L = 0), unfavorable (R L > 1), linear (R L = 1), and favorable (0 < R L < 1). The obtained R L values in this study were in the range of 0−1, thus confirming that the Tesla'Ssorb is a suitable adsorbent for the adsorption of S 8 under defined testing conditions. Graphical presentations for Langmuir and Freundlich equilibrium isotherms of S 8 adsorption onto Tesla'Ssorb (linear and non-linear models) are given in Figures 9-11. Correlation coefficients (R 2 values) were used to compare the suitability of used isotherm models. Since the obtained R 2 values for Langmuir isotherm model were higher in comparison with Freundlich isotherm (Table 2), it can be concluded that the Langmuir isotherm fits better. The experimental data, in both linearized and non-linearized fitting curves, were fitted better to the Langmuir isotherm model, as shown in Figures 9-11. Furthermore, it can be observed that maximum adsorption capacity (q m ) increases with increasing temperature, due to the endothermic nature of adsorption. This was further confirmed by the analysis of thermodynamic parameters.

Adsorption Thermodynamics
The thermodynamic parameters of elemental sulfur removal from mineral insulating oils using Tesla'Ssorb provide information about: the adsorption process (whether it is spontaneous or not), the nature of interactions between the Tesla'Ssorb and S 8 molecules, and the conditions for reaction which will provide the highest rate of S 8 removal. The adsorption mechanism can be determined though thermodynamic parameters, such as change in Gibb's free energy (∆G • ), enthalpy of adsorption (∆H • ), and entropy (∆S • ). The change in the Gibbs free energy (∆G • ) was calculated from the distribution coefficient (Table 3). The relationship is shown in Equation (8), as follows: where ∆G is the calculated change in the Gibbs free energy, R is the universal gas constant (8.314 J·mol −1 ·K −1 ), and T is the absolute temperature in Kelvin, respectively.  K d is the distribution coefficient which can be calculated from the concentration of S 8 adsorbed at equilibrium (q e ) and S 8 concentration in oil at equilibrium (C e ), as follows: The relationship of (∆G) to enthalpy change (∆H) and entropy change (∆S) of adsorption as well as the relationship between (∆G) and the (lnK d ) are shown in Equations (10) and (11), respectively: Substituting Equation (8) into Equation (10) gives the Equation (11): The values of ∆H and ∆S are determined from the slope and the intercept of the linear plot of lnK d vs. 1/T respectively. The plots of lnK d vs. 1/T, calculated at three temperatures, 328, 338, and 353 K, are shown in Figure 12.

Adsorption Kinetics
Pseudo-first-order, pseudo-second-order, and Weber-Morris intraparticle diffusion models were used to analyze the adsorption kinetics of S8 on Tesla'Ssorb.
A linear form of pseudo-first model is: where k1 is the pseudo-first order kinetics rate constant (min −1 ), calculated from the linear plots of log(qe − qt) vs. t. The negative values of Gibbs free energy for S 8 adsorption indicated that the process was spontaneous, and more spontaneous and thermodynamically favorable at higher temperatures, for various initial concentrations of S 8 .
The calculated values shown in Table 3 suggest an adsorption based on combined physical adsorption and chemical reaction, i.e., chemisorption process [34].
The positive ∆H values determined for all initial S 8 concentrations have indicated that the reaction is endothermic in nature. The randomness at the solid-liquid boundary layer during the adsorption can be described by standard entropy changes, ∆S. The ∆S values were positive for all initial concentrations, which indicated an increase in randomness of the solid-liquid phases at the boundary layer and a high affinity and probability of reaction between the S 8 in mineral oil and metal ions on Tesla'Ssorb.

Adsorption Kinetics
Pseudo-first-order, pseudo-second-order, and Weber-Morris intraparticle diffusion models were used to analyze the adsorption kinetics of S 8 on Tesla'Ssorb.
A linear form of pseudo-first model is: where k 1 is the pseudo-first order kinetics rate constant (min −1 ), calculated from the linear plots of log(q e − q t ) vs. t. Pseudo-second order has a linear form as follows: where k 2 is the pseudo-second order kinetics rate constant (g mg −1 min −1 ). The Weber-Morris intraparticle diffusion model investigates the intraparticle uptake of adsorbate and the pore diffusion in adsorption. The linear form of this model is: where k id (mg/g min 1/2 ) is the intraparticle diffusion adsorption rate constant, and C is a parameter related to the boundary layer's thickness. Calculated kinetic parameters summarized in Table 4 for different initial S 8 concentrations at different temperatures indicated that the data fitted to the pseudo-second-order model provided a higher value of correlation coefficient R 2 (from 0.9972 to 0.9999) in comparison to the pseudo-first-order and intra-particle diffusion kinetic models. Moreover, the calculated equilibrium adsorption capacity (q e,cal ) from pseudo-second-order kinetics is closer to the experimental value (q e,exp ). This is also in line with calculated values of maximum adsorption capacities of Tesla'Ssorb (q m ) at different temperatures using the Langmuir isotherm (Table 2). This implies that the adsorption may be governed by chemical adsorption. The graphical presentations for pseudo-first, pseudo-second order kinetic (linear and non-linear), and intra-particle diffusion models at three different temperatures are given in Figures 13-15  The graphical presentations for pseudo-first, pseudo-second order kinetic (linear and non-linear), and intra-particle diffusion models at three different temperatures are given in Figures 13-15 respectively.     The operating conditions, such as temperature and initial concentration of S 8 , have the influence of the value of k 2 . This decreases as the initial S 8 concentration increases, since it takes a longer to reach equilibrium (Figure 14). At higher concentrations, lower adsorption rates are obtained, due to high competition for the Tesla'Ssorb surface active sites.
The rate controlling steps that occurred during the adsorption process were determined using the Weber-Morris model. Figure 15 shows that the plot was not linear over the whole time range and had two linear segments, indicating that intra-particle diffusion was not the only rate-controlling step and that the sorption process is more complex and involves several diffusion resistances.
At all temperatures, the slope of the first portion was higher than the slope of the second portion, indicating that the S 8 molecules diffused more quickly through the boundary layer than through the pores of the adsorbent. Initially, the S 8 molecules were rapidly adsorbed on the surface of the Tesla'Ssorb due to chemical interactions between S 8 and incorporated Ag ions on active sites. Later on, after saturation was reached, the S 8 molecules diffused into the internal surfaces of the adsorbent particles, further followed by the establishment of equilibrium. The operating conditions, such as temperature and initial concentration of S8, have the influence of the value of k2. This decreases as the initial S8 concentration increases, since it takes a longer to reach equilibrium (Figure 14). At higher concentrations, lower adsorption rates are obtained, due to high competition for the Tesla'Ssorb surface active sites.
The rate controlling steps that occurred during the adsorption process were determined using the Weber-Morris model. Figure 15 shows that the plot was not linear over the whole time range and had two linear segments, indicating that intra-particle diffusion was not the only rate-controlling step and that the sorption process is more complex and involves several diffusion resistances.

Activation Energy (E a )
The activation energy (E a ) determines the correlation and dependence of reaction between S 8 and silver with temperature. In absorption processes, E a is the minimum amount of energy needed for interaction between the adsorbate, i.e., elemental sulfur and incorporated silver ions in Tesla'Ssorb.
The activation energy, E a (kJ mol −1 ) for S 8 adsorption was determined from the Arrhenius equation, at three different temperatures (328, 338 and 353) K. The Arrhenius plot was developed by plotting the lnk 2 vs. 1/T (in Kelvin). The Arrhenius equation's linearized form is provided in Equation (15): where k 2 is the pseudo-second-order rate constant (g·mg −1 ·min −1 ), E a is the activation energy of adsorption (kJ mol −1 ), and A is the Arrhenius constant. The value of E a defines the type of adsorption, which may be physical or chemical. The activation energy for physical adsorption is usually no higher than 4.2 kJ mol −1 , due to the weak forces that occur in physical adsorption. Chemical adsorption is more specific and involves much stronger forces than in physical adsorption. There are two types of chemical adsorption, activated and non-activated as not frequently present. According to a finite activation energy in the Arrhenius equation (8.4-83.7 kJ mol −1 ), activated chemical adsorption occurs when the rate changes with temperature. The rapid occurrence of a non-activated chemisorption suggests that the activation energy is close to zero [35]. A linear plot between lnk 2 and 1/T shows a straight line ( Figure 16). The activation energy (E a ) for the removal of the highest concentration of S 8 (153.1 mg/kg) was found to be 27.1 kJmol −1 , suggesting that S 8 adsorption by Tesla'Ssorb is controlled by an activated chemical adsorption.
The value of Ea defines the type of adsorption, which may be physical or chemical. The activation energy for physical adsorption is usually no higher than 4.2 kJ mol −1 , due to the weak forces that occur in physical adsorption. Chemical adsorption is more specific and involves much stronger forces than in physical adsorption. There are two types of chemical adsorption, activated and non-activated as not frequently present. According to a finite activation energy in the Arrhenius equation (8.4-83.7 kJ mol −1 ), activated chemical adsorption occurs when the rate changes with temperature. The rapid occurrence of a non-activated chemisorption suggests that the activation energy is close to zero [35]. A linear plot between lnk2 and 1/T shows a straight line ( Figure 16). The activation energy (Ea) for the removal of the highest concentration of S8 (153.1 mg/kg) was found to be 27.1 kJmol −1 , suggesting that S8 adsorption by Tesla'Ssorb is controlled by an activated chemical adsorption.  The obtained result of activation energy is in good correlation with thermodynamic parameters, thereby confirming the postulated mechanism of chemical reaction (Reaction 1) and high affinity, i.e., reactivity of molecules S 8 and silver.

Conclusions
The adsorption performance of new mesoporous material with spherical clusters of silver nanoparticles, HSA-Tesla'Ssorb, was investigated in this study. The Tesla'Ssorb was synthesized, characterized, and tested for the removal of S 8 in wide range of concentrations (from 79.7 to 153.1 mg/kg). The removal of sulfur from the mineral oil is obtained by chemical reaction of Ag ions (deposited on the Tesla'Ssorb surface), with sulfur, followed by neutralization of acidic by-products with ammonia. FESEM, EDS, and XRD analyses confirmed the presence of incorporated silver in Tesla'Ssorb, which is essential in the removal of S 8 from the mineral oils.
The results of adsorption experiments on a pilot scale unit, in terms of the effect of contact time and equilibrium absorption capacity of Tesla'Ssorb, match very well with the results previously reported. This confirms that the pilot scale unit is a very good tool for predicting an adsorbent's on-site performance, even if a high concentration of S 8 is present in the oil (153.1 mg/kg). The effect of temperature was pronounced at the highest initial concentrations of S 8 , confirming that 353 K can be used as an optimal temperature for the efficient removal of S 8 from the oil, in a wide range of concentrations.
Experimental results agree very well with the Langmuir adsorption isotherm model. The thermodynamic parameters from the study indicated that the adsorption of the S 8 is spontaneous and endothermic.
According to the kinetics studies, it was observed that the adsorption of S 8 is very fast in the beginning and decreases towards approaching the equilibrium. Kinetics data were best fitted by the pseudo-second order model, indicating that the rate-limiting step in the adsorption process is chemisorption.
The results of the intra particle diffusion kinetic model indicated that the adsorption process could be divided into two stages. The fast adsorption in the first stage was driven by boundary layer diffusion while the slow adsorption in the second stage was attributed to intraparticle diffusion where the equilibrium adsorption was achieved.
The activation energy (E a ) calculated from Arrhenius equation also indicated activated chemical adsorption of S 8 onto Tesla'Ssorb. Compared to various conventional adsorbents used in oil reclamation processes, a low amount (up to 3 wt.%) of Tesla'Ssorb was verified to have high capacity to remove variable concentrations of S 8 from mineral insulating oils in a short treatment time.